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Abstract 



We study the finite temperature symmetry behaviour of O(Ni) x 0(N2) scalar models 
on the lattice and we prove that at sufficiently high temperatures and in arbitrary di- 
mensions their full symmetry is always restored or, equivalently, that the phenomenon of 
Symmetry Non Restoration which, according to lowest order perturbation theory, takes 
place in the continuum version of these models, does not occur on the lattice. 
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The high temperature behaviour of relativistic field theories has been the subject of an intense research 
since the early works by Kirzhnitz and Linde Weinberg g] and Dolan and Jackiw ||. The results 
of these investigations, show that in a "typical" case the symmetry of the vacuum increases when the 
temperature is raised and thus, in spontaneously broken theories, internal symmetries are gradually 
restored when the system is heated up. Thus, in the context of Grand Unified Theories and Cosmology 
the fact that the universe might have been in phases with different symmetry properties at different 
stages of its evolution, undergoing a series of phase transitions in the process of cooling down, would have 
significant consequences, like for instance the creation of topological defects via the Kibble mechanism 

Nevertheless, as noticed by Weinberg ||, O(Ni) x 0{N2) theories may present an "atypical" symmetry 
behaviour, where by atypical, we mean that either the symmetry is not restored at high temperatures or 
that an exact symmetry of the low temperature theory becomes broken at higher temperatures. These 
phenomena, known as Symmetry Non Restoration (SNR) and Inverse Symmetry Breaking (ISN) are 
in fact two aspects of the same problem, the only difference between them being whether or not the 
symmetry is broken in the zero temperature theory. To lowest order in perturbation theory, the existence 
of SNR or ISB is related to the possibility of having negative Debye masses. This can be achieved 
in multi-scalar theories, provided that some of the fourth-order couplings are taken negative and large 
enough in absolute value (but also small enough as to produce a bounded potential). As on the other 
hand, the scalar sector of most extensions of the Standard Model and Grand Unified Theories is rather 
undetermined, it turns out that SNR and ISB are not so atypical as one would first suppose. 

Indeed, it has been recognised that the ideas of SNR and ISB can have interesting phenomenological 
implications. Very recently for instance, the phenomenon of SNR has been used by Dvali, Melfo and 
Senjanovic j5| to suggest that the monopole problem might not exist in some GUT's, by arguing that the 
monopole-producing phase transition might have never occurred (for an earlier implementation of this 
idea see [||). Also in connection to monopoles, ISB is the basis of the proposal by Langacker and Pi 0, 
which states that a period of broken U(l) em can cure the monopole problem. Other no less interesting 
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applications of SNR and ISB concern the breaking of the CP symmetry ]|], the domain wall problem 
H , Baryogenesis , Inflation 1 1 and P, Strong CP, and Peccei-Quinn Symmetries . ISB and SNR 
have also been considered in |l3[ ]. 

Appealing as the idea of SNR and ISB may be, no consensus has yet been reached on whether they 
really correspond to a true physical effect or rather to an artifact of ( may be, lowest order?) perturbation 
theory. This second point of view stems from the fact that when non-perturbative approximations are 
used to study the symmetry behaviour of these models flifl , it is found that symmetry is invariably 
restored at high temperature. Moreover, even staying within the realms of perturbation theory, it has 
been shown that the inclusion of next to leading order effects in the calculation of thermal masses 
tends to reduce the region of parameter space in which SNR and ISB occur. Recently the gap 
equations used in p"5| have been rederived by a more detailed analysis based in the Cornwall- Jackiw- 
Tamboulis effective potential [jlTj at finete temperature Jl8| . An independent study, which also encodes 
some non perturbative information through the Effective Action Technique [ ft9| , has been carried out 
[2C|] . For the small values of the scalar self-couplings which were considered in J2(| , it was found that the 
corrections to the lowest order perturbative computation are small. While this is not in disagreement 
with the results of |]l5| , which also predict small corrections when these couplings are small, we think that 
an extension of the analysis of [20[| to include large couplings may be unavoidable, since, as pointed out 
in realistic models may require large scalar self-couplings, due to the presence of gauge interactions 
which conspire against SNR and ISB. 

Due to the conflicting results which emerge from the perturbative, semi perturbative and non per- 
turbative methods we mentioned above, we think that Lattice Field Theory might result to be a useful 
setting to study SNR and ISB. Some work along this lines has been already done in j2l) where it is shown 
that, when an approximation based on the constraint effective potential p2[ is made, symmetry is always 
restored in less than four dimensions. 

In this letter, we will study SNR and ISB on the Lattice without making approximations and in 
arbitrary dimensions and we will prove that symmetry is always restored, at sufficiently high temperatures, 



2 



for 0(N%) x O(A^) models. In doing this, we will closely follow the steps of a theorem by King and Yaffe 
on symmetry restoration for O(N) models on the Lattice. It is worth mentioning that while King and 
Yaffe's result for O(N) models confirms the result obtained from perturbation theory, the result we 
present here exactly contradicts the lowest order perturbative calculation. 

The continuum model we initially consider is a global 0(N±) x 0(N2) scalar theory in (d+ 1) euclidean 
dimensions described by the action, 




where 4>i = (<Pi > ' ' ' ■> 91 ) is an iVj-component real scalar and | <pi | 2 = X}j=iW ) 2 - The condition of 
boundedness for the potential constrains the coupling constants to satisfy the relations: 

A, > , 4AiA 2 > A 2 . (2) 

This condition allows for positive values of A and, as it immediately follows from the one-loop computation 
(|], if A is such that, 

A > !±p„ (3, 

the 0(N\) x 0(N2) symmetry is necessarily broken to 0(N\) x 0(N2 — 1) at high T. The symmetry at 
lower temperatures depends instead on the signs and magnitudes of the masses m 2 . 

We will now prove that when the model (|l|) is defined on a discrete lattice of points, the full O(Ni) x 
0(iV 2 ) symmetry is restored at sufficiently high temperatures, for all values of the parameters in the 
action. 

We consider an anisotropic hypercubic (d+ l)-dimensional lattice AeTxE. Here, £ is an infinite d- 
dimcnsional hypercubic lattice accounting for space, while T is a finite one-dimensional lattice consisting 
of iV T points, accounting for the finite euclidean time axis. We assume a priori distinct spacings Ax and 
At for S and T respectively. Thus: 

A ee T x E = {x = (xq, x) : xq = hqAt x = (niAx, • • • , n^Ax) ; 

no = 1,- • -,iV T ; rii e Z i = 1,- • -,d} . (4) 
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The temperature of the system is then given by: 



T = 



1 



iV T Ar 



(5) 



The lattice version of the model is described by the action: 

1 ( 1 



S = ^A-r(A:r) rf ]T 



xeA 



2 \ (At) 2 



d 1 \ 
(x + e T Ar) - ^(as) | 2 + V -r-r^ | <f>i(x + e k Ax) - &(z) | 2 



1 



+ -m 2 | ft (a) | 2 +Ai(| ft(x) | 2 ) 2 



A | 4>i{x) | 2 | 2 (x) | s 



(0) 



where e T and are unit vectors pointing along the time and fc-th spatial directions respectively and the 
fields 4>i(x) satisfy periodic boundary conditions in the time direction: 



4> t (i + AtN t , x) = fa (xo , x) 



(7) 



It is convenient to introduce the dimcnsionlcss quantities: 



4>i(x) = (Ax) 2 4>i(x), a 



At 



mj = (Ax) 2 m 2 , Aj = (Ax) 3 ~ d \ l7 A = (Ax) 3 - d A. 
In terms of them the action (||) reads: 

5=X){Z] \ \ \ \<l>i{x + e T AT)-(t>i{x) | 2 +53 a l <t>i(x + e k Ax) - <t>i(x) \ 2 \ 



+ iam 2 | &(x) | 2 +aA l (| ft (a;) | 2 ) 5 



aA|^(x) | 2 |</» 2 (x) | ; 



(8) 



(9) 



It is also convenient to measure the temperature in units of the inverse of the lattice spacing Ax. We 
thus define the "lattice temperature" 



T L = (Ax)T = 



1 1 



A way to test the symmetry of the system is to consider the two-point functions: 
(Uz)-^(w))=Z- 1 / f II II d ^ x n JM-hi™) exp(-S) , 



(10) 
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where Z is the partition function and z and w are two spacelike separated points of A. In the broken 
phases, one or both of the correlators ( |TT| ) have a non-vanishing limit M? for infinite separations: 



Mf = lim (4>j(z) ■ 4>j(w)) ^ for some j (12) 

\z — w\ — >oo 

while in the symmetric phase both these limits vanish. Now, we will prove that, for all values of the mass 
parameters and coupling constants in ([)]), and in any dimension d, there exists a temperature T£ above 
which the correlators ( fill) decay exponentially with the separation | z — w \ : 



\< c ^eM-M{T L )\z-w\] T L >T* L . (13) 
In this equation c is a numerical constant, M(Tl) is a function such that: 

M(T L ) > for T L > T* L ; lim M{T L ) = oo (14) 



and v is a constant independent on Tj,. Once (g_3j) is proved, it then will follow that above T£ both limits 
Mj vanish and the full O(Ni) x 0(N 2 ) symmetry is restored. 
Before proceeding further, a few comments are in order: 

a) the temperature T£ provides only an upper bound for the true critical temperature T£ of the 
symmetry- restoring phase transition. As pointed out already in P3| , this bound is not expected to have 
the correct dependence on the bare coupling constants, for large values of the latter. Consequently, based 
on this bound, no conclusions can be drawn for the critical temperature in the continuum limit, whenever 
this limit exists. 

b) the bound we are going to derive, like the one in [^3) for the O(N) models, depends on a only 
through the temperature Tl, — aN T . This will allow us to take the continuum limit in the time direction 
a — * 0, aN T — > const in a straightforward way. 

We now turn to the proof, which is, as we mentioned before, a simple extension of King and Yaffe 
theorem on symmetry restoration for O(N) models p3| . The intuitive idea behind it is that, at high 
temperature, the system can be thought of as a collection of oscillators sitting on the sites of the spatial 
lattice S, such that no order is possible. Guided by this idea, we will write the action, the partition 
function and the correlators in a way which will turn to be useful to derive the bound ( ^3|) 



We will start by rewriting the action (0) as: 



s = E s * - E v ^ - 



(15) 



where E* is the set of links in S and 



x £T i 



while 



-z- | (f>i(xo + At, x) - (f>i(xo,x) \ 2 +~a(mf + 2d) \ 4>i(xo,x) | 2 +aA»(| 4>i(xo,x) | 2 ) 2 
2a 2 



aA | </>i(x ,x) \ 2 \ (t) 2 (x , 

x £T 



(16) 



F w = ^E E«i^o,i) + toj)i 2 , (17) 

i £T i=l,2 

where x and y are the end points of the link I. Now, the first term on the r.h.s. of ( |l5| ) precisely describes 
a set of uncoupled oscillators located at the sites of S, while the second sum provides an interaction 
among them. Notice also that the V(l) is positive, a property which will be important in what follows. 
Upon defining the measure 

dfl = Yl dfi x 

where 

dfi s = z~ 1 Y[ Yl d(j>i(xQ,x) exp[-5 s ] . (18) 

x eT i=l,2 

and z is a coefficient that normalises dfi s to one, the partition function and the correlators can be written 



as 



Z = I dfi exp 



Eno 



(19) 



■j(z) ■ 4>j(w)) = Z 1 / dfj,(f)j(z) ■ 4>j(w) exp 



If we now define: 



expV(0 = l + p(0 



(20) 



(21) 



and 



tf(Y) = [df^Yl P(l) 
J leY 



(22) 



then 

z= fdn n [i+ P {i)]= Yi K( y) ■ ( 2a ) 

J ze£* res* 

In the same way, defining 

Ki (Y)= fdn J~(z) ■ (w) J] p(l) , (24) 

J l BY 

we then have: 

{J-{z) ■ ^{w)) ^ Z- 1 ]T K 3 {Y) . (25) 

FCS* 

Here and in ( p3| ) the sum is over all subsets Y of £*. 

Let now W be the connected component of Y that contains w and let X = Y — W. As the action 5 5 
is invariant under (pi — > —<f>i, the only non- vanishing terms in the sum ( p^ ) are those for which z is in W. 
Together with the fact that the measures dfi s are normalised to one, this allows us to write: 

(^(i) •&(«>)>= E E_| K ( X ) ( 26 ) 

W connected 

where W denotes the closure of W, namely the set of links in S* that share an end-point with some link 
in W. 

Up to here, we have only rewritten Z and the correlation functions in a convenient way. We shall now 
start to look for bounds for these quantities. The positivity of V(l) comes now into play in a crucial way, 
as it immediately implies: 

* J_ W K{X) - w^m^vw s 1 w c • (27) 



which in turn implies the bound: 



(M*) E K AW). (28) 

w c s* 

connected 

The next step is then to find a bound for Kj(W). First of all, we observe that: 

| P (l) |=| expV(0 - 1 |<| V(l) | exp V(0 . (29) 



Moreover, by Schwarz's inequality, we have 



V ^ ^EEH &(so,2) I 2 +a | ^i(xoj) I 2 ] , 

x eT i=l,2 



(30) 



which implies the other bound 



q(x) 



Uv(i)< £ n £ £ 

ieW {g(S)}S£S lx GTi=l,2 



(31) 



Here the sum is over all possible choices of the non-negative integers q(x) < 2d which vanish for all x £ W 
and are such that 

£<?(*) =1^1 (32) 
where | W | denotes the number of links in W. Use of the inequality: 



finally implies the bound: 



^Wl< £ / ^(l^-WI 2 + I^HI 2 ) 



(33) 



n 

2G£ 



(9(5)} 

E E a i ^(^0)*) 

i eT t=i,2 



q(x) 



exp 



p ( x ) E E a I ^( x o,a;) 

i eT i=l,2 



(34) 



The factor p(x) in the above formula represents the number of links in W that have x as endpoint. 
Obviously 

p{x) < 2d, Vx . (35) 

The important feature of eq.(|34|) is that Kj(W) is now bounded by a sum of products of independent 
one-dimensional integrals. The next move is to bound the latter by Gaussian integrals. To this purpose, 
we define: 



dvx = z 1 1 JJ JJ d<p l (x ,x)exp <-- 

x £T i=l,2 [ 



E E ~ I ^*( x o + Ar,x) - ^i(x ,a;) | 2 +a/i 2 | <j>i(x ,x) | 2 
— L — a 

x eT 1=1,2 

(36) 

In this equation, z' is a normalisation factor, while /j, 2 represents a variational parameter, whose value 
will be fixed at the end such as to provide the best bound. For simplicity, we have introduced a common 



parameter /i 2 for both fields 4>i, while in principle a better bound could be obtained by letting distinct 
values. The derivation of the bound would in this case be sightly more involved but the final bound 
would not be qualitatively different from the one obtained with one parameter only. 

Let us now go back to ([mJ): each term in the sum is, as we said, a product of independent one- 
dimensional integrals, I Sl one for each x € W. There are now two cases to be considered: whether x 
is distinct from z and w or not. In the first case, with the aid of (|36|), we can write the corresponding 
one-dimensional integral I s as: 

N x 

j- = - 

X „_ 



with 



N x = / dv x \ E E a I 4>i{ x o,x) I 2 expT x [p(x)] 

\x £Ti=l,2 J 

D x = J dv x exp T x [p = 0] 



(37) 
(38) 



and 



t*w)\ = E E 

x eT i=l,2 



i/i 2 - irn 2 +p{x) - 2d) a | ^(xq^x) | 2 -oA»(| 4>i(x ,x) | 2 ) 2 



x er 



(39) 



Now, it is trivial to check that T x [p(x)} is bounded by 



Tx\p(x)\ < 



4AiA 2 - A 2 



A , ( i^ 2 - l -m\ + p(x) - 2d\ + Ai Q/i 2 - im 2 + p(x) - 2d\ + 



A, 



(40) 



where we have used Ylx eT a — Pl- We then find that: 



N x < exp(A) / E E a I M x o,x) | 2 

J \x £Ti=l,2 J 



q(x) 



(41) 



To proceed further, we will need to use the following bounds for the moments of a Gaussian distribu- 



tion: 



dug | Mx l o\x) | 2 ••■ I </h(x$",x) | a < (2p- 1)!! 



(p) ^ |2< 



Pl vm 2 



(42) 
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du s | <f>i(x ,x) | 2 > 



fay 2 



dv S | tf>i(x ,x) | 2 | fa(x ,x) | 2 > o2 4 

PlA 4 



to write: 



7V S < exp(yl) ^ ^ ^ y ' dv s ^2 a\ 0i(x o ,a;) | 2 ^j ^ ^ a | 2 (^o,a;) I 



q(x)—r 



q{x) 



<(4d-l)!!^ 



r=0 



r 
g(x) 



(4d- 1)!! 



v, l^ + A 



(N 1+ N 2 ) ( 4 + — 



1 , h 



g(x) 



exp„4 



As for D 2 , we can use Jensen's inequality (sec p4j) to bound it as: 



exp.4 



(43) 



(44) 



D 2 >exp/d^[0]> 



exp 



E 

i=l,2 



^ / 1 .2 1 2 



M 2 V 2 



A* - o m i - 2d ) - 3Ai/?L 



~fa^ 



exp(-B) 



Using (45) together with (44) we obtain 
J* < (4d-l)!! 



(N 1+ N 2 ) ( 4 + — 
A* A* 



q{x) 



x cxp(.4 + B) 



(45) 



If instead x coincides with either zorui the previous bound has to be multiplied by a factor of 

PL VA- 4 A* / 

Since now the only sites cc for which p(x) or g(x) are different from zero are those which are endpoints 
of some link in W and since there are less than 2 | W | such sites, we can use the previous bounds and 
eg. (^)to say 



Kj(W) |< 2(4d+l)- 



Nj ( 1 , fa 



fa VA* 2 A 1 



E 

{9(2)1 



M A 4 



|W| 



[(4d-l)!!] 2 l w/ lexp(2 | | (.4+B)) 

(46) 



The number of terms in the sum above is less than (2d) ' ' (see |23| for a proof) and thus we have: 



c(iVi + N 2 ) ( ^ + — 
\ l A* 



exp(2|W|(„4 + B)), (47) 
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where c = 4ci 2 [(4d — l)!!] 2 . If we now take: 



fa ' 



(48) 



where v is independent on /3l, there will be a fa such that, for (3l < (3l 



Hfa<l fJ 2 >\ m} | i = 1,2 



(49) 



These relations, together with (B5J) then imply: 



K 3 (W) |<4(4d+l)^=L= 



2c(N l +N 2 )\l — 
v 



\w\ 



x exp { 2 | W 



2d{N x + N 2 )J & + /\+ A2+ ^ + -(^XxNl + 12A 2 7V 2 ) 
V v 4AiA2 — A- v 



Ifajs 



exp \-\W\ M'(fa) 



(50) 



It is clear that lim^^^o M'(fa) = oo and thus there exists a (5* L < fa such that, for 0l < fa, M'((3l) > 0. 
The steps to derive the final bound ( |T3| ) from eq. (^0|) are exactly the same as in ]23| and we omit repeating 
the short proof here. (/3£) _1 thus represent our bound for the true critical temperature of the symmetry- 
restoring phase transition: 



n< fa 



(51) 



fa depends of course on v, d, Ai, A and mf and one might further exploit the freedom in the choice of v 
in order to get the best bound; as the analogue best bound for f3* L in the O(N) case [|3j does not have 
the right behaviour for large values of the bare couplings and as we believe that the same will occur in 
our case, we will not explicitly show it here. 

The above proof can be extended to the case when the O(Ni) x 0{N2) symmetry is gauged (partly 
or completely), along the lines of |2J|, and we refer the reader to that paper for the details. 

We have proved that, in arbitrary dimensions and at sufficiently high temperatures, the full 0(N\) x 
0(N%) symmetry is restored , or equivalently, that the phenomenon of Symmetry Non Restoration does 
not occur on the lattice. Notice, though, that the implications of our proof regarding Inverse Symmetry 
Breaking are weaker, in the sense that we can not exclude the possibility of this phenomenon taking 
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place at intermediate temperatures. We can only state that if an ISB phase transition occurs at a given 
temperature, then the system will necessarily undergo a symmetry restoring phase transition at higher 
temperatures. 

Whether these results are relevant for the continuum, when the continuum limit exists, depends on 
the behaviour of the critical temperature of the symmetry-restoring phase transition, when the lattice 
spacing Ax is taken to zero. Analogously to the O(N) case, our bound for T c diverges when the bare 
couplings become large and thus it is not possible to say based on this bound, if T c remains finite in the 
continuum limit. We would like to stress, once more, that this behaviour of the bound is not peculiar 
to the case examined here, but appears also in the O(N) case in p3| , for which one knows that T c has 
a finite continuum limit. An answer to this question using Monte Carlo simulations is under current 
investigation. 
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